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EXPLOSIVE BREAKUP 




• Explosive breakup of jettisoned spacecraft segments during entry is a relatively 
poorly understood phenomenon 

• ATV Jules Verne Breakup 

http: / /www.youtube.com/watch?v=8os300hZLfE 

• Ejected fragments may collide with main entry vehicle and cause catastrophic 
damage that can result in mission failure 

• Highly probabilistic event dependent upon jettisoned segment geometry, entry 
trajectory, component makeup, etc. 

- Difficult to predict breakup altitude or direction of ejected fragments 

- Makes deterministic assessment very challenging 
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PREVIOUS APPROACH TO DEBRIS RE-CONTACT 




• Typically the re-contact problem has been approached using large set Monte 
Carlo runs that examine closest approaches to the entry vehicle 

• However, re-contact events are rare, therefore even a large data set without 
conjunctions does not provide strong evidence for improbability 

• Computationally expensive to generate very large sets of Monte Carlo data, 
especially when adding dispersions 


• Hinders the ability to rapidly optimize separation sequences (jettison angles, 
separation burn magnitude, etc.) 
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MESH-BASED APPROACH 



• Sought to develop a method that would aid in optimizing separation 
trajectories by minimizing the risk of debris re-contact to an entry vehicle in 
the case of an explosive breakup of a jettisoned segment 

• Required a statistical approach to determine the likelihood of a debris strike 
given a certain spatial distribution of debris and separation trajectory 

• Non-Homogeneous Poisson Process 

- Debris strike probability proportional to the cumulative debris flux along 
a particular trajectory 

• Allows more robust analysis than simply considering the closest approach 
data from several thousand Monte Carlo runs 

• Chosen because it is a relatively simple probabilistic model to explain and 
offers a reasonable approximation of the risk of re-contact 
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SIMULATION OF FRAGMENTS 




• Explosive breakup simulated by imparting an ejection velocity to each fragment 
in a random direction while assuming: 

- Instantaneous breakup 

- Constant drag model for debris pieces (no lift) 

- No collisions between fragments 

- No debris demise 2 xl0 "' 


• Empirical evidence from observed 
breakups can be used to determine 
fragment ejection velocities and 
ballistic coefficients 

- NASA Standard Breakup Model 



• T+30s After Breakup 

• T+5s After Breakup 
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NASA STANDARD BREAKUP MODEL 



• Uses empirical data from gathered 
from observed orbital breakups 

• Gives the number, area to mass 
(A/M) ratio and corresponding 
delta-v for fragments from the 
explosive breakup of spacecraft 
and upper stage rocket bodies. 

• Ballistic coefficient easily 
calculated from A/M ratio, 
assuming a C D = 1 


Delta- V vs. Ballistic Coefficient for Explosive Breakup 
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DEBRIS DISTRIBUTION MAPS 




• Region of space around the jettisoned 
body divided up into a mesh grid 
describing volumetric cells 


Debris Density Map T+3s after breakup 
(Hove spherical and Cartesian equivalent as well) 


• Spatial debris density computed 
using the weighted number of 
fragments inside each cell at any 
point in time 
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Spatial Debris Density, Pieces/m 




Altitude, 




DEBRIS DISTRIBUTION MAPS 


Evolution of Debris Density Map with Time (w/ debris fragments overlaid) 
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DEBRIS FLUX 


• As the entry vehicle moves along its trajectory relative to the jettisoned 
body and its fragments, it moves through the meshed region of space 


• Use weighted spatial density inside a cell to calculate debris flux, (p D at 


each point in time along a particular trajectory 
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DEBRIS DENSITY ALONG A TRAJECTORY 



Debris Density/Relative Velocity Along a Particular Generic Separation Trajectory 



o' ^ — 1 'o 

195 200 205 210 215 220 

Elapsed Time, sec 


2011 AIAA Modeling & Simulation Technologies Conference 


10 


08/09/2011 


Debris-Vehicle Average Relative Velocity, m/s 






APPLICATION OF STATISTICAL MODEL 



• Using this debris flux we can determine the intensity parameter A, as a 
function of time along a separation trajectory 

• A(t) is the product of the debris flux and the cross-sectional area of the 
entry vehicle, taken as a constant equal to its maximum circular area, A v 

Mt) = <p D (t)-A v 

• This can then be integrated over the interval [t v t 2 ] to give the cumulative 
intensity parameter, A and therefore, the total probability of re-contact 
along a particular trajectory 

1 2 

A = jA(r)dr 

a 

P{X> 1} = 1 - ^ A 
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ASSUMPTIONS OF STATISTICAL MODEL 




• A debris strike is a rare event such that, 

P{X > 1} « 0 

• The probability of a debris strike at any time is an independent random variable, 

P{X 1 \X 2 } = P{X 1 } 

• Then frequency with which a strike is expected to occur is dependent only upon 
the debris flux ((p D ) in a region of space or volumetric cell 

• Debris fragments are small relative to the entry vehicle 

• The entry vehicle can be modeled as a sphere with uniform cross-sectional area 
in all directions 
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Cumulative Re-contact Probability 




CUMULATIVE RE-CONTACT PROBABILITY 



Cumulative Re-contact Probabilities Along Four Separation Trajectories 
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• Useful information from an operations perspective, to predict when an 
entry vehicle is 'under threat' from debris 
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SELECTION OF MESH CELL SIZE 



Adaptive method for selection of cell size allows cells to grow with the debris cloud 


• Reduces discretization errors 
owing to the sparsity of 
fragments, particularly at 
greater distances from the blast 

•More effective at capturing 
the evolution of the fragments 
with time 

- Volume of the cells grow as 
the uncertainties in 
fragment position grow 


Adaptive vs. Fixed Mesh Cell Size 
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ADVANTAGES OF THIS APPROACH 



• Intuitive explanation of the likelihood of a debris strike - 'good' and 'bad' 
regions of space are easily visualized 

• Computationally efficient for optimizing entry vehicle separation trajectories 

- Only requires one set of debris data to run several trajectory cases 

• Robust enough to handle a wide variety of separation maneuvers 

- Posigrade, retrograde, out-of-plane burns 

• Allows capturing of additional data including relative velocity and ballistic 
coefficients of debris pieces 

- Capability to then use this to additionally calculate a consequence of 
impact metric, and hence the risk of a debris strike 
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QUESTIONS 
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BACK-UP SLIDES 
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SELECTION OF MESH TYPE 



• Mesh type can lead to discretization error where debris 
cloud is not well described by the mesh 

- Boundaries of cells don't 
match the debris cloud 

- Insufficient fragments 
modeled for the size of the 
cells 
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Total Re-contact Probability for Three Mesh 
Types as Separation Burn Time is Varied 


• Cartesian mesh tends to 
exaggerate these errors due to 
fixed volume cells 

• Careful selection of size and 
shape must be made to avoid 
these errors in the model 
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m DERIVATION OF POISSON PROCESS PROBABILITY 



The probability of k events occurring can be modeled using a Non-Homogeneous 
Poisson Process according to the following formula: 


P{X=k } 



-A 


where, h 

A = JA(r)dT 

t 2 

If we consider the probability of k = 0 (no strikes), then the probability of one or 
more strikes occurring is given by: 


P{X=0} = 



-A 


=>P{X =0} = e~ A 


=> P{X > 1} = l — e~ A 
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